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ABSTRACT 

PTF09dav is a peculiar subluminous type la supernova (SN) discovered by the Palomar Tran- 
sient Factory (PTF). Spectroscopically, it appears superficially similar to the class of subluminous 
SN1991bg-like SNe, but it has several unusual features which make it stand out from this population. 
Rs peak luminosity is fainter than any previously discovered SN1991bg-like SN la (Mb ~ —15.5), but 
without the unusually red optical colors expected if the faint luminosity were due to extinction. The 
photospheric optical spectra have very unusual strong lines of Sc II and Mg I, with possible Sr II, 
together with stronger than average Ti II and low velocities of ~ 6000 km s^ 1 . The host galaxy of 
PTF09dav is ambiguous. The SN lies either on the extreme outskirts (~ 41 kpc) of a spiral galaxy, 
or in an very faint (Mr > —12.8) dwarf galaxy, unlike other 1991bg-like SNe which are invariably 
associated with massive, old stellar populations. PTF09dav is also an outlier on the light-curve- 
width-luminosity and color-luminosity relations derived for other sub-luminous SNe la. The inferred 
56 Ni mass is small (0.019 ± 0.003M Q ), as is the estimated ejecta mass of O.36M . Taken together, 
these properties make PTF09dav a remarkable event. We discuss various physical models that could 
explain PTF 09dav. Helium shell detonation or deflagration on the surface of a CO white-dwarf can 
explain some of the features of PTF09dav, including the presence of Sc and the low photospheric 
velocities, but the observed Si and Mg are not predicted to be very abundant in these models. We 
conclude that no single model is currently capable of explaining all of the observed signatures of 
PTF09dav. 

Subject headings: supernovae: general — supernovae: individual (PTF09dav) 



1. INTRODUCTION 

Type la supernovae (SNe la), thermonuclear explo- 
sions of accreting carbon-oxygen white dwarf stars, form 
a fairly uniform and homogeneous class of events. This 
has inspired their application as cosmological standard- 
izable candles, tracing the cosmic expansion history to 
high redshift. However, the class of "sub-luminous" 
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SN la events show quite different photometric and 
spectroscopic properties than those mainstream events 
used for cosmological studies, with potentially differ- 
ent progenitor models. Though the discovery of these 
rare, faint SNe la has historically been challenging, the 
archetypal sub-luminous SN 1991bg (jFilippenko et al.l 
1992; ILeibund gut et al.l I1993D has been complemented 
with the discove ry of many s i milar events (e.g., see 
comp i lations of iBranch et al.l 119931: iGarnavich et al 
| 2004t iTaubenberger et al.l l2008t iGonzalez-Gaitan et al 
l2010f h Modern SN searches are revealing ever more 
extreme examples of faint SNe of all types, the na- 
ture of which , in some cases, remains uncertain (e.g 
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The key characteristics of the 1991bg-likc sub-luminous 
SN la population are a faint peak absolute magnitude 
(~ 2 magnitudes fainter than a normal SN la) with 
a red optical color, a light curve that rises and falls 
~ 40% faster than a normal SN la, a strong trend of 
being lo cated in old, E/SO galaxies w ith large stellar 
masses (jHowelll l200lt iNeill et al.l I2009D . and 'cool' op- 
tical maxim um light spectra dominated by strong Ti II 
absorption (Filipp enko et al.l 119921: IMazzali et"aTTll997l: 
IBranch et al.|[2006D with some what lower expansi on ve- 
locities than in normal SNe la dBenetti et al.ll2005f). The 
class of SNe la similar to SN2002cx (|Li et all 120031 ) ap- 
pear to form a further, separate sub-class of faint events, 
with a host galaxy distribution favoring later morpho- 
logical types comp ared to that of SN1991bg-like events 
(|Folev et all 120091: IValenti et al.ll2009D . Some events in 
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the sub-class, such as SN2008ha, ar e fainter still than 
1991bg-like SNe (|Folev et all [20091: IMcClelland eTall 
2010), with even lower ejecta velocities (and kinetic en- 
ergies). Recent observations also demonstrate the exis- 
tence of further faint or fast types of SN events, where 
the formal SN type classification may provide mis-leading 
information about the physical nature of the explosion 
(jPoznanski et al.ll20"lfl iPe'rets et al.l 120 lObllal ). 

Various physical models and scenarios have been pro- 
posed to explain sub-luminous SN1991bg-like events. 
Whether they form a distinct physical group from normal 
SNe la, with different progenitors and explosion models, 
or whether they lie at the extreme end of a continuous 
distribution, but with lower 56 Ni masses and hence tem- 
peratures, is unclear. No conclusive evidence for any 
one progenitor model has so far been found. The na- 
ture of the progenitors of SN2002cx-like events is even 
less clear, with possible models ranging from the direct 
coll apse of a massive s t ar to a black hole (a "fallback" 
SN; iValenti et all [20091: iMoriva et aL|[201Qh. to the pure 
deflagration of a w hite dwarf (e.g. iBr anch et all [200l 
IPhillins et afll200l . 

In this paper we describe the supernova PTF09dav. 
This SN was dis covered as p art of the Palomar Transient 
Factory fPT^ Ilau et al.|[20ql lLaw et al.ll2009h . a five 
year project surveying the optically transient sky. PTF 
is designed to sample a large fraction of the optical tran- 
sient population, including both new types of events, as 
well as statistically complete samples of known transient 
types. In particular, the high-cadence, survey depth, and 
wide-area make PTF i deal for the study o f fast, sublu- 
minous SN events (e.g. lKasliwal et al.ll2010T ). PTF09dav 
is one such event - a very subluminous SN la with pe- 
culiar (and so far unique) spectral properties. In §[2] we 
introduce PTF09dav together with its basic photomet- 
ric, spectral, and host galaxy data. We analyze the light 
curve in § [31 deriving estimates of its peak luminosity 
and related parameters. § 0] details our spectral analy- 
sis, including the identification of the principal elements. 
We discuss various physical interpretations of our data 
in § El In a companion paper, we discuss the late-time 
light curve and nebular spectroscopic data for this event 
(Kasliwal et al., in prep.). Where relevant we assume 
iJ =70 km s" 1 Mpc" 1 . 

2. OBSERVATIONS 

PTF09dav was discovered as a transient event at 
R ~ 19.5 by the Palomar Transient Factory (PTF) on 
2009 August 11.3 (all dates in this paper are UTC), lo- 
cated at a = 22 h 46 m 55!15, S = +21°37'34'.'l (J2000; 
Fig. [TJ . The event was isolated with no apparent host 
galaxy in the PTF reference images. A spectrum was 
taken with the William Herschel Telescope (WHT) us- 
ing the Intermediate dispersion Spectrograph and Imag- 
ing System (ISIS) on 2009 August 14.0 as part of the 
PTF SN la key project. The R316R (red arm) and 
R300B (blue arm) gratings were used, together with the 
5300 dichroic, giving a wavelength coverage of ^3400- 
8100A. Spectral comparisons were performed at the tele- 
scope using the SU PERFIT spectral matching code (e.g., 
iHowell et al.l 12005), and the spectrum was initially rec- 
ognized as being superficially similar to the sublumi- 



nous SN la SN1991bg (|Fili D penko et al.lll99l . with deep 
Ti II, Si II, and O I absorption, located at a redshift of 
z ~ 0.04 (Fig. [2]). Given the faint discovery magnitude 
and relatively low redshift, PTF 09dav was recognized as 
a potentially rare and interesting sub-luminous SN event, 
and was triggered for further study, which we summarize 
in the remainder of this section. 

2.1. Spectra 

Following the initial WHT classification spectrum, two 
further spectra were taken during the photospheric phase 
of the SN. The first was taken with the Palomar Hale 
200-in t elescope (P200) plus Double Beam Spectrograph 
fDBSP: [Qke fc Gunnl Il98l on 2009 August 20.5. The 
6001/4000A (blue) and 1581/7500A (red) gratings, and 
D55 dichroic, were used, giving a wavelength coverage of 
~3500-9500A. The second spectrum was taken using the 
Keck-I telescope and the Low-Re solution Imaging Spec- 
trometer (LRIS; lOke et al.l [1991 on 2009 August 25.5 
- a wavelength coverage of -3400-10000A was achieved 
using the 4001/3400 A grism (blue), 4001/8500A grating 
(red), and the 560 dichroic. A final spectrum taken 
during the nebular phase of the SN was taken on 2009 
November 11.4, and will be presented in a companion 
paper (Kasliwal et al. in prep.). 

The three photospheric spectra were reduced using the 
same custom-written pipeline based on standard proce- 
dures in I RAF and I PL broadly followin g the reduction 
procedures outlined in lEllis et al.l (|2008t ). including flux 
calibration and telluric feature removal. All spectra were 
observed at the parallactic angle in photometric condi- 
tions. "Error" spectra are derived from a knowledge 
of the CCD properties and Poisson statistics, and are 
tracked throughout the reduction procedure. The spec- 
tra are also rebinned (in a weighted way) to a common 
dispersion, and joined across the dichroic. The three 
spectra are shown in Figs. [2] and [3] 

2.2. Host galaxy 

The host galaxy of PTF09dav is somewhat ambigu- 
ous (Fig. PTF09dav lies 56.8" (projected) SE of 
a potential host galaxy (r ~ 17.6; a = 22 h 46 m 52?9, 
5 = +21°38'21'.'6, J2000). A Lick spectrum taken on 
2009 October 25th shows this galaxy to have strong neb- 
ular emission lines (Ha, Hf3, O III, O II) and be located 
at a heliocentric redshift Zh i = 0.0371 ± 0.0002 (well 
within the Hubble flow), consistent with the estimated 
SN redshift. At the position of the SN, z^ c \ = 0.0371 cor- 
responds to z — 0.0359 in the CMB rest-frame (z cm b), or 
a distance modulus of 35.99 ± 0.06 (the error assumes a 
host galaxy peculiar velocity uncertainty of 300 km s -1 ). 
This would give a projected physical separation between 
PTF 09dav and the center of this galaxy of 40.6 kpc. 

There is no evidence for a host galaxy at the positio n of 
the SN. Using 129 "DeepSky"E3 f N ugent et al.ll2009tt im- 
ages taken with the Palomar-QUEST survey over of the 
field of PTF 09dav, we have constructed a deep RG610 
image (a longpass r + i + z filter). Using a 4" diam- 
eter aperture at the position of the SN, we detect no 
host galaxy with a 3-cr limiting magnitude of m = 23.2 
in the RG610, or M = -12.8 at z = 0.0359. Such a 
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faint absolute magnitude would be lower than any de- 
tected host galaxy of a core collapse SN found by PTF 
ijArcavi et al.l [2010) . Deeper Keck imaging at the posi- 
tion of PTF 09dav is underway (Kasliwal et al., in prep.). 

2.3. Light curve data 

Photometric monitoring in gri filters com menced using 
the r obotic Palomar 60in (P60) telescope (jCenko et al.l 
120061 ) on 2009 August 15.3 and on the robotic Faulkes 
Telescope North (FTN) on 2009 August 18.4, comple- 
menting the rolling PTF search on the Samuel Oschin 
48in telescope (P48) in i?-band. The resulting multi- 
color light curve of PTF 09dav is shown in Fig. 0] As 
PTF 09dav has no detected host galaxy at the position of 
the SN, no subtraction of a reference image is required, 
simplifying the light curve measurement. We measure 
the SN photometry using a point-spread-function (PSF) 
fitting method. In each image frame, the PSF is deter- 
mined from nearby field stars, and this average PSF is 
then fit at the position of the SN event weighting each 
pixel according to Poisson statistics, yielding a SN flux 
and flux error. 

We perform a flux calibration to the Sloan Digital Sky 
Surv ey (SPSS) photomet ric system, close to the AB sys- 
tem (|Oke & Gunnlll983f ). using P48 and P60 observa- 
tions of SDSS fields made on the same night as obser- 
vations of PTF09dav. We use an average calibration 
from photometric nights in August and September 2009, 
which provides calibrated tertiary standards in the field 
of PTF09dav. This calibration is then transferred to 
other epochs, and to the FTN observations, using aper- 
ture photometry of these tertiary standards to determine 
an relative flux-scaling factor for each epoch, typically 
measured to <1%. In the P60/FTN g, r and i filters 
the color terms to the SDSS system are very small. The 
P48 color term is larger, and we include color, extinc- 
tion, and color-extinction terms. The r.m.s. of the color 
term fits is 0.02-0.03 mag (with a color term in r — i of 
~ 0.22 mag), although we present magnitudes in the nat- 
ural P48 system and do not apply this color term to our 
photometry. Finally, we fit the same PSF used for the 
SN flux measurement to the calibrating field stars, and 
apply a correction to go from the aperture in which the 
flux calibration is determined, to the PSF fit in which 
the SN is measured. 

Although we could perform our entire analysis using 
magnitudes defined in the AB system, rest-frame SN la 
properties in the literature are typically given in the Vega 
system. For ease of comparison to previous work, we 
convert our calibrated light curves into this system when 
performing the light curve analysis, and all quoted mag- 
nitudes are given in the Vega system. The offsets from 
Vega (m ve ga) to AB (tuab) magnitudes are -0.12, 0.13, 
0.19, and 0.35 for the g, r, R, and i filters respectively 

(?Tivega 

= toab - offset). Our light curve data can be 
found in Table [IJ presented in counts rather than mag- 
nitudes to preserve information on non-detections. 

3. LIGHT CURVE ANALYSIS 
3.1. Method 

Estimating the peak luminosity of SNe in a given 
bandpass requires an interpolation between observed 
data points at the time of maximum light, followed 



by a fc-correction back to some standard rest-frame fil- 
ter of interest. The interpolation step is best per- 
formed using a smooth light curve template in the ob- 
served filter which is believed to represent the photo- 
metric evolution of the SN; the fc-correction requires 
a time series spectral energy distri bution (SED; e.g. 
iNueent et all 120021: iHsiao et al.ll2007t) . To parametrize 
the light curve of PTF 09dav, we fit a time series 
SED to the obs erved photometry using the SiFTO 
light curve fitter (jConlev et al.l I2008Q developed for SN 
la cosmology studies. We replace the standard SN 
la spectral template^ with on e based on the sub- 
luminous SN la events S N1991bg dFilippenko et al.lll992l: 
Leibundeut et al.lll993ft and SN1999bv (|Garnavich et~atl 
2004). Ideally one would use a SED template which ex- 
actly matches the spectral class of the object under study. 
For PTF 09dav, the first object of its type, no such tem- 
plate exists, so we use our sub-luminous template which 
has many common spectral features and similar broad- 
band colors. 

SiFTO works in flux space, manipulating a model of 
the SED and synthesizing an observer-frame light curve 
from a given spectral time-series in a set of filters at 
a given redshift, allowing an arbitrary normalization in 
each observed filter (i.e., the absolute colors of the tem- 
plate being fit are not important and do not influence the 
fit). The time-axis of the template is adjusted by a di- 
mensionless relative "stretch" (s) factor (where the input 
template is defined to have s = 1) to fit the data. Once 
the observer-frame SiFTO fit is complete, it can be used 
to estimate rest-frame magnitudes in any given set of 
filters, provided there is equivalent observer-frame filter 
coverage, and at any epoch. This is performed by ad- 
justing the template SED at the required epoch to have 
the correct observed colors from the SIFTO fit, correct- 
ing for extinction along the line of sight in the Milky 
Way, de-redshifting, and integrating the resultant SED 
through the required filters. This process is essentially a 
cross-filter k-correction, with the advantage that all the 
observed data contribute to the SED shape used. As an 
illustration, the size of the k-correction for a SN1991bg- 
like SN la to go from observer g (r) at z = 0.0371 to 
rest-frame B (V) is ~ —0.16 mag (~ —0.24 mag) at max- 
imum light, and ~ —0.28 mag (~ —0.48 mag) at 10 days 
after maximum. 

The function used to adjust the template SED can 
either be an interpolating spline (in which case a per- 
fect match in color is attained) or, if the reason for the 
physical difference in the colors between the SN and the 
template is u nderstood, a functio n such as a standard 
dust law (e.g. iCardelli et al.lll989f) or other wavelength- 
dependent color law (e.g. iGuv et all 120071 ). The latter 
case is more useful when extrapolating beyond the wave- 
length range of the observed filters, although it gives a 
poorer match between observed and fit colors. 

3.2. Fit results 

Our SiFTO light curve fit is in Fig. g] The fit is good - 
the reduced \ 2 of the fit is 1.02 for 40 degrees of freedom 
(i.e. x 2 = 40.8/40) - using a standard SN la template 
in SiFTO gives a substantially worse reduced \ 2 °f 4.3. 

16 http : //supernova. lbl .gov/$\sim$nugent/nugent_templates .html 
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Maximum light in the rest-frame -B-band was on 2009- 
08-08.9 ± 0.3 days (or an MJD of 55051.9), i.e. the SN 
was discovered some 3 days after maximum light. Inter- 
estingly, although the fit to the data around maximum 
light is excellent, at phases beyond ~+30 days there ap- 
pears to be significant excess i-band flux in the SN com- 
pared to the simple prediction of the template with the 
i-band light curve appearing to plateau, while the decay 
in the other filters continues to follow the sub-luminous 
template. This could be due to emerging nebular emis- 
sion from the [Ca II] doublet at A7291 and A7323, with 
stronger lines at that phase than in our SN1991bg-like 
sub-luminous template, and more similar to lat er time 
spectra of calcium-rich SNe (jPerets et al.ll2010"bT ). 

As PTF 09dav has no detected host galaxy at the SN 
position, we can cross-check our photometric calibration 
using photometry synthesized from our WHT spectrum. 
At the time of the WHT spectrum (+4.9 days), SiFTO 
gives a light curve color of g — R = 1.17 ± 0.05, com- 
pared to a spectral color from the WHT spectrum of 
g — R = 1.26. Given the uncertainties in flux-calibrating 
long-slit spectra (e.g. differential slit losses which will 
redden the SED), as well as in the P48 R and P60 g 
filter responses used for synthetic photometry, this is an 
excellent level of agreement suggesting an accurate pho- 
tometric flux calibration. A similar agreement is found 
for the later P200 spectrum. 

Our reported magnitudes are in t he Vega s ystem , tak- 
ing our BVR filter responses from Besscll (1990), and 
correcting for Milky Way extinction assu ming Ry = 3.1 
and a color e xcess of E(B — V~) mw = Q,Q44 (Schlcgc l et al.l 
119981 ) with a iCardelli et all (j 19891) extinction law. We 
measure peak rest-frame absolute magnitudes at the time 
of B-band maximum light of M B = -15.44 ± 0.05, 
M v = -16.01 ± 0.05 and Mr = -16.26 ± 0.05 (er- 
rors are statistical uncertainties propagated through the 
light curve fit only) using a CCM law in SiFTO to 
match the template to the observed fit colors; we adopt 
these numbers throughout this paper. For compari- 
son, using a spline in SIFTO in place of the CCM law 
gives M B = -15.54 ± 0.05, M v = -15.98 ± 0.09 and 
M R = -16.24 ±0.05. 

One potential concern in the light curve fit is that the 
P60/FTN g coverage, corresponding to rest-frame B, did 
not start until ~ 6 days after maximum light - and so 
the -B-band (and y-band) maximum light estimate is an 
extrapolation (based on the light curves synthesized from 
the SED template) rather than an interpolation as with 
the P48 i?-band. Our key assumption is that the light 
curve of PTF 09dav is well-represented by a "stretched" 
version of our subluminous spectral template. While this 
is clearly true in the P48 R data where the fit before and 
after maximum light is robust, in g the data are not 
sufficient to test this assumption. 

3.3. Comparison to other events 

A comparison of PTF09dav to other SN1991bg-like 
subluminous SN la events can be found in Table [2] 
and Fig. [6l With no commonly accepted definition 
for what constitutes a SN1991bg-like event, we fit a 
combi nation of eve n ts ma tching the cool classification 
from iBranchetafl (120091) . events with AM 15 (B) > 
1.8 fr om iTaubenberger et alj (|2008f) and (|Hicken et all 
2009), and other individual events from the literature 



(e.g. IKasliwal et all l200l . We fit all the public SN 
la photometry with SiFTO and the same subluminous 
template to ensure a consistent comparison. All magni- 
tudes are corrected for Milky Way extinction, and we 
have assumed distance moduli either calculated from 
Zcmb when z cm b > 0.01, or taken from a redshift- 
independent estimate (together with an uncertainty) 
from the NASA/IPAC Extragalactic Database (NEEEZI) 
in other cases. We have also propagated through a 
300 kms -1 peculiar velocity error for those SNe la in 
the Hubble flow. We only consider SNe with adequate 
photometry for a reliable light curve fit, and those with 
Mb > —18 at maximum light, in our analysis. 

The SNe show expected trends between stretch, color 
and M b (Fig[6l) which are consistent with previous 
studies dGarnavich et al1l2004l ; ITaubenberger et aLll2008t 
iGonzalez-Gaitan et alTl2010f ): Fainter SNe are redder, 
and with smaller values of stretch. (Note that the 
stretches here are measured relative to our subluminous 
template and are not on the same system as stretches 
of normal SNe la.) Clearly PTF09dav is an outlier in 
these relations, being considerably fainter for its color 
and stretch than other subluminous events: PTF09dav 
is ~ 3.5 mag under-luminous compared to normal SNe la 
(Mb — —19.1), ~ 1.5 mag fainter than a typical sub- 
luminous SN la, and around ~ 0.5 mag fainter than 
SN2007ax, the faint est subluminous SN la identified 
(jKasliwal et al.ll2008f ). 

No correction for host extinction (if any is present) has 
been applied to the magnitudes of PTF09dav (or the 
other SNe in Table [2]), and thus host extinction will im- 
pact Fig. E For PTF 09dav, (B - V) msK = 0.57 ± 0.08, 
and when using the CCM law in SiFTO to match the 
subluminous template to the fit colors for PTF 09dav, 
a small additional E(B - V) of 0.08 ± 0.04 is required. 
Though it may be tempting to interpret this at face value 
as a reddening for the SN (and derive similar correc- 
tions for the other SNe in Table [2]), such an approach is 
likely to be dangerous. Generally the color-luminosity 
variation in SNe la is poorly understood. There is con- 
siderable evidence for an intrinsic color-luminosity rela- 
tion which is currently impossible to disentangle from 
the effects of dust unless the SN is particularly well ob- 
served. Thus interpreting a red color as evidence of ex- 
tinction may not be correct. A second consideration is 
that our subluminous template may not be appropriate 
for zero extinction, as it is based on real SN photome- 
try Therefore we simply note that the broad-band colors 
of PTF09dav are similar to other sub-luminous SN la 
events and show no evidence for significant host galaxy 
extinction. 

We also estimate the value of the iPhillipsI (|1993| ) de- 
cline rate parameter Ami5(-B) (the amount in magni- 
tudes in the rest-frame -B-band that the SN declines in 
the 15 days following maximum light) as Amis(-B) = 
1.87 ± 0.06. The error is a combination of propagat- 
ing the error in the stretch parameter, combined with 
the uncertainty in peak B-band brightness. Numer- 
ous studies have fit the relation between Awi5(-B) and 
Mb for subluminous SNe la (e.g. JGarnavich et al.ll2004l : 
ITaubenberger et alJ l2008f). A ccording to the relation 
of ITaubenberger et all (|2008h . our Am 15 (B) = 1.87 

17 http : //nedwww. ipac . caltech.edu/ 
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should correspond to a peak B absolute magnitude of 
— 17.44, nearly two magnitudes brighter than we observe 
{Mb — —15.5 corresponds to Am±5(B) ~ 2.15 according 
to the Taubenberger relation). This is consistent with 
PTF09dav as an outlier in Fig. [BJ However, Am 15 (_B), 
unlike stretch, can be misleading in these comparisons. 
In the fainter examples of subluminous SNe la, the tran- 
sition from the fast initial decline to the slower late-time 
decline o ccurs before —15 day s, which distorts the rela- 
tion (e.g. fKasliwal et al.ll2008l ). 

3.4. Rise-time and energetics 

The subluminous template we used in the SiFTO fit 
has a rise-time (r r ) from explosion of 13.0 days to -B-band 
maximum light, and 14.0 days to bolometric maximum 
light. PTF09dav has a light curve "stretch" factor of 
0.86 ± 0.02 compared to this template, giving a rise-time 
to bolometric peak of 12.0 ± 0.3 days (statistical error 
only). We also estimate the rise time by fitting the P48 
R data in the early rise-time region using the analyti- 
cal equation f(t) — a(r + r r ) 2 (Ries s et al.lll999[ ). Here, 
r = t/(l + z), where t is the observed time, and a is a 
normalizing constant. We find rise-times of 12.0 ± 2.9 
when fitting the first two epochs (including the first non- 
detection), and 14.5±1.1 when including the next epoch. 
These are consistent with our estimates from the tem- 
plate light curve fit. The latter measure includes data 
from near- maximum light, where the assumption of the 
simple rise-time model may no longer hold, and should 
be considered a strong upper limit. 

Using our light curve fits, we can also estimate the 
bolometric luminosity, Lboi, of PTF09dav. Estimating 
the bolometric luminosity is never straightforward as full 
ultraviolet through infrared (or UVOIR) light curves are 
rarely observed. Instead we measure the luminosity in 
a given band together with the colors from our light 
curve fits, and use the template to derive a bolomet- 
ric luminosity estimate (effectively a bolometric correc- 
tion). Such techniques work well for SNe la events near 
maximum light as the peak of the luminosity output is 
in the optical. Again, as with the light curve fit, we 
use the template spectrum to calculate our bolometric 
light curve. In principle we could use the observed spec- 
tra, but the lack of an observation at maximum light, 
and the lack of A > 7700 data in the spectrum nearest 
maximum, mak es this difficu l t. Ou r technique broadly 
matches that of iHowell et a l. (2009). We take the tem- 
plate spectrum on the day of maximum light, adjust it 
using a CCM law to have the correct color, normalize 
to the peak B-band absolute magnitude, and integrate 
the resulting spectrum. The peak bolometric luminosity 
(occurring ~ 1 day after i3-band peak) is estimated to 
be L^ ol = 5.6 ± 0.4 x 10 41 crgs- 1 (Table [2]). The im- 
plicit bolometric correction, defined as BC = M^ol — My 
where M(, j is the bolometric magnitude, is ~ 0.2 mag, 
comparable to other sub-luminous events SN la events 
(jContardo et al.|[2000h . 

The ejecta mass M c j and kinetic energy E'km of 
PTF 09dav can be estimated from its rise time and ejecta 
velocity v e j using M e j ex v e jT^ and £kin oc v^t^ (e.g., 
lArnettl Il98l iFolev et all [2501 ) . Comparing PTF 09dav 
to a "typical" SN la with a kinetic energy of 10 51 ergs, a 
rise-time of 17.4 days (jHavden et al.|[2010l) . and a photo- 



spheric velocity of 11,000 kins" 1 (|Benetti et al.l 12005ft . 
and assuming the opacities of both are the same, we 
estimate an ejecta mass of 0.36M Q and kinetic energy 
of 8 x 10 49 ergs. (Here, we have used an ejecta veloc- 
ity of 6000 km s -1 for PTF09dav estimated from the 
spectral fitting in § 31) These are low values; a typ- 
ical ejecta mass for a SN la ranges from ~ 0.5Af Q 
for a SN1991bg-like event, up to th e Chandrasekhar 
mass (1.4M(p) for normal SNe la (e.g. iStritzinger et al.l 
120061 : iMazzali et al.l [20071 ) . A similar low ejecta mass 
was derived for th e calcium-rich SN lb SN2005E 
(jPerets et al.l l2010bT) . and even lower ejecta masses 
have been e stimated for faint examples of SN2002cx - 
like events (IFolev et al.l l2009t iMcClelland et all 12010ft 
and for SN2002bi dPoznanski et al.l l2010ft and SN2010X 
(jKasliwal et al.ll2010D . 

If the light curve is predominantly powered by the de- 
cay of 56 Ni, the 56 Ni mass, Mm, can be e stimated from 
a knowledge of the rise-time r r and Lboi (|Arnett][l982). 
At maximum light, Mm is given by MNi=iboi/aS'(r r ) 
(Arnett's rule), with S the radioactive luminosity per 
solar mass of 56 Ni from the decay to 56 Co to 56 Fe eval- 
uated at maximum light, and a the ratio of bolometric 
to radioac tive luminosit i es, ne ar unity. We adopt S a s 
defined in IHowell etaD (|2009f ) (see also iBranchl 11992ft . 
which includes the effect of r r , and take r r and its error 
from our light curve fits with a conservative additional 2 
day uncertainty in the rise time error estimate to match 
the range found above. In practice, a is likely to de- 
viate from unity. For normal SNe la, the light curve 
may peak at a luminosity that exceeds the instantaneous 
radioactive luminosit y (a > 1) due to a falling tempera- 
ture and opacity fe.g. lBranchlll992l ). On the other hand, 
for smaller ejecta masses some gamma rays may escape 
the ejecta without b eing thermalized (i.e., a < 1). Us- 
ing the approach of Kasliwal et al. (201(J) we estimate 
the optical depth of PTF09dav at maximum light to 
be ~ 70 indicating an effective trapping of emitted 7- 
rays. For simplicity we therefore take a = 1, giving 
Mffi — 0.019 ± 0.003M Q . Under the same assumptions, 
SN 2008ha produced 0.0029M Q of 56 Ni, and SN 2007ax 
0.038 ± O.OO8M0 - see Table [2] for estimated M Ni for all 
the 1991bg-likc subluminous SNe. In these calculations, 
we estimated r r for each event by multiplying the stretch 
by the rise-time of the subluminous template (14 days). 
We emphasize that this calculation assumes the entire 
light curve to be powered by the decay of 56 Ni. 

In summary of this section, the Mb, Lboi, Mnj and 
ejecta mass M e j of PTF 09dav are all unusually low. 
PTF 09dav is one of the faintest subluminous SNe la yet 
discovered, and, while not as faint as SN2008ha, it is the 
faintest of the SN1991bg-like sub-class. 

4. SPECTRAL ANALYSIS 

We now turn to the interpretation of the PTF 09dav 
spectra. The spectra of PTF 09dav show many curious 
features compared to other su b-luminous events (F ig. [2]) . 
The comparison to SN2005H (|Taubenberger et al.l 2008) 
is particularly instructive, as that spectrum was taken at 
an identical light curve phase (+5d). Despite the resem- 
blance to sub-luminous SNe la, with obvious Ti II, Si II, 
and O I absorption, other features in the spectra do not 
appear in other sub-luminous events (Fig. [2]). In particu- 



G 



Sullivan et al. 



lar, absorption features at 3960A, 5400A, 5540A, 5720A, 
6480A, 6700A, and 6890A do not, at first glance, im- 
mediately correspond to features in other sub-luminous 
SNe la. We attempt to identify these lines using an im- 
plementation of the SYNDW code, and then discuss the 
results. 

4.1. Spectral fitting 

When no detailed ah initio SN explosion model 
provides an immediate explanation for observations, 
parametrized spectrum synthesis provides the first step 
towards their construction. The idea is to use simplified 
radiative transfer calculations to directly fit SN spectra. 
A good fit constrains explosion models through inter- 
pretive spectral feature identification, with the main re- 
sult being the detection or exclusion of specific chemical 
elements. The velocity distribution of detected species 
within the ejecta can also be constrained. 

In our analysis we make use of SYNAPPS (Thomas et 
al., submitted). The physical assumptions SYNAPPS uses 
match those of the well-known SYNOW code (|Fisherll2061l . 
so findings are restricted to identification of features 
and not quantitative abundances. But where SYNDW is 
completely interactive, SYNAPPS is automated. This re- 
lieves the user from tedious, iterative adjustment of a 
large number of parameters (over 50 variables) to gain 
fit agreement, and assures more exhaustive searching of 
the parameter space. SYNAPPS can be thought of as the 
hybridization of a SYNDW-like calculation with a parallel 
optimization framework, where spectral fit quality serves 
as the objective function to optimize. 

4.2. Spectral Line Identifications 

We run SYNAPPS on our three photosphcric spectra, 
and the resulting fits are shown in Fig.[5j We identify the 
following lines common to subluminous SN la spectra: 
O I, Ca II, Si II, Ti II, Fe II, Co II. The Ti II lines 
are particularly strong relative to other SN1991bg-like 
events. We also identify Sc II, Na I, and Mg I, as well 
as evidence for Sr II and possibly Cr II. There was no 
evidence for S II or Mg II - both degraded the quality of 
the fits, S II around 5500A, and Mg II in the red. 

Though unusual, the presence of Sc II seems robust. As 
well as the two features at 5400A and 5550A (caused by 
A5527 and A5658), the 6490A feature (caused by A6604), 
together with improved fits between 4000-5000A, pro- 
vides additional confirmation. S II has strong lines in 
this region, but cannot be responsible for the two fea- 
tures at 5400A and 5550A (observer frame) - the ratio 
of the two line wavelengths does not match S II, and the 
velocity would be inconsistent with the other elements 
(Fig. [TJ . S II is common in 1991bg-like events (although 
the lines typically become weaker and disappear after 
maximum light), and its non-detection here is therefore 
surprising, but it does have a higher ionization energy 
than the other elements, which may point to lower tem- 
peratures. 

Although Sc II has been includ ed in the SYNDW 
fits of faint SN2002cx-like SNe la (|Folev et all [200l 
IMcClelland eFldl l2010h . it has never been robustly de- 
tected in SNe la. Sc II has been observed (together 
with the other s-process elements Ba II and Sr II) in 
SN1987A (|Williamslll987T; iMazzali eTaI1 IT992). and also 



durin g the plateau pha s e of some lower-lumino sity SNe 
IIP (jFassia et al.l ll~998t iPastorello et~aT1 l200l . Mg I 
(with no indication of Mg II) is also surprising. As 
both Sc II and Mg I are only seen at low temperatures 
(jHatano et al.lll999Q . even lower than Ti II, this provides 
further evidence for the low temperatures of PTF 09dav. 

Na I A5892 is also highly probable (Fig. [7]), and is typ- 
ically seen in subluminous events after maximum light, 
although it is particularly strong in PTF 09dav - Fig. [3j 
This feature could also be He I A5876 line, and we cannot 
unambiguously confirm Na I over He I (Fig. [7|). Includ- 
ing He I in the fits does better match the 6900A feature 
(produced by A7065), but it introduces other features 
in the SYNAPPS spectra which are not seen in the data, 
such as 6500A (caused by A6678). However, SYNAPPS 
does not explicitly account for non-local thermodynamic 
equilibrium effects whic h are known to be important for 
modeling He I features (|Lucvl[l99"l] ) - LTE excitation in 
SYNAPPS can overestimate the strength of A6678 to A5876 
and A7065 {B ranc 1: 2003). Thus, all we can say is that at 
least one of, and possibly both, Na I and He I are present 
in the spectra of PTF 09dav. 

There is also reasonable evidence for the presence of 
Sr II (via absorption near 3960A) , with weaker evidence 
for Cr II (primarily through absorption at 4770 A), all 
shown in Fig. [7] Sr II is identified through a blended 
single line (A4078 and A4216), and has no other strong 
optical features. As with Sc II, Sr II is also only seen 
at low temperatures (jHatano et alJ 119991 ) . Cr II is of- 
ten weakly present in SN1991bg -like SNe (Mazzal i et al.1 
Il997t iTaubenberger et al.ll2008| ). Ba II, an s-process ele- 
ment often seen with Sr II, was not needed in the fits. 

Given the apparent low temperatures, we also tested 
the inclusion of S I and Si I, but neither was required in 
the fits (their abundances came out to zero in SYNAPPS). 
Indeed, S I and Si I would add lines not seen in our spec- 
tra if their abundance were non-zero. We also experi- 
mented with V II and C II - there was no evidence for 
either in the fits, although they also did not degrade the 
fit results by adding extra lines not seen in the spectra. 

4.3. Spectral measurements 

We measure ejecta velocities from the blueshift of the 
absorption minima of the P-Cygni line profiles. This 
process can carry a number of uncertainties, particularly 
the blending of neighboring lines, so we use the SYNAPPS 
results which effectively average over several lines. The 
line velocities evolve from ~ 6100 kms -1 (+4.9 days) to 
~ 5100 kms" 1 (+11.2 days) to ~ 4600 kms" 1 (+16.0 
days), and the major lines all share a common velocity. 

The strength of the Na I (or He I) line increases with 
time, and our tentative Sr II detection also becomes 
stronger. Ti II, already stronger than in most SN1991bg- 
like events, also increases in strength with phase. By 
contrast, the evidence for Fe II, as well as Cr II, has dis- 
appeared by +16 days; those lines are not needed in the 
later-epoch SYNAPPS fits 

We also estimate the Nugent ct al. (|1995f ) Si II ratio, 
7\L(Si), defined as the ratio of the depth of the Si II fea- 
tures at A5972 and A6355. We calculate this ratio as 0.35, 
smaller th en the values typical of su b-luminous SNe la 
of ~ 0.65 (jTaubenberger et all [20081. and naively indi- 
cating a higher temperature than in those events. How- 
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ever, the presence of strong emission from the Na I (or 
He I) line blue-wards of the A5972 feature will signifi- 
cantly affect this measurement, decreasing the depth of 
the A5972 feature and hence artificially decreasing the 
observed ratio. We also measure the (pseudo) equiv- 
alent widths of these two Si II features, which can 
be used to differentiate sub-types of SNe la usin g the 
W(6100)/W(5750) plane (IBranch et al.ll2006L 12009^ . Us- 
ing the notation of IBranch et all (|2006D .~Wf6100) (the 
A6355 line) is 45A and W(5750) (the A5972 line) is lOA. 
These measurements are both significantly lower than in 
SN1991bg-like events (typical values of ~100A and ~45A 
respectively) , and therefore lie well outside the region in 
W(6100)/W(5750) space co ntaining other cool SNe with 
significant Ti II absorption (Bra nch et al.|[2009l ). 

5. DISCUSSION 

PTF09dav is an unusual supernova. Though super- 
ficially similar to a subluminous SN la, it is consider- 
ably fainter at Mb — —15.5 (without being redder), and 
is one of the fastest SNe la observed. The peculiarity 
extends to its spectra, which have low ejecta velocities 
(~ 6000 kms" 1 ) and unusual spectral lines, including an 
unambiguous detection of Sc II - which we believe to 
be unique to PTF09dav among SNe la. The origin of 
Sc is not clear. Though it is the decay product of 44 Ti, 
the half-life is ~ 60 years, which makes that an unlikely 
source for the PTF09dav maximum light spectra. 45 Ti, 
with a half-life of ~ 3 hours, or higher and more unstable 
Ti isotopes, are more plausible if they are synthesized in 
large enough quantities. Some 45 Ti is predicted to be 
synthesize d in the delayed d etonation SN la explosion 
models of Ma eda et al.l ()2010D . but only with very small 
mass fractions. The low temperatures also probably con- 
tribute to the enhanced visibility of Sc. If real, the Sr II 
detection is also unusual, and presumably originates from 
the s-process in the progenitor star. 

PTF 09dav does share some similarities with SN2005E 
and the sub-class of calcium-rich Type Ib/c SNe 
(jPerets et alJl2010b[ ). SN2005E was faint (M B = -14.8), 
as are all calcium-rich SNe, with a similar low ejecta mass 
to PTF09dav (but higher ejecta velocities). SN2005E 
had a Ca rich spectrum with no prominent intermediate- 
mass C/O-burning products (e.g., Mg, S, Si), and al- 
though PTF 09dav is also Ca/Ti-rich, unlike SN2005E it 
does show the presence of Si and Mg lines. Perhaps more 
interestingly, the environments are very similar. The en- 
vironment of PTF 09dav is discrepant compared to other 
subluminous SNe la: PTF09dav is either lying on the 
extreme outskirts of a spiral galaxy (~ 41kpc from the 
center, projected distance), in contrast to the massive 
E/SO galaxies that form the bulk of the SN1991bg-like 
host population, or is located in a very faint, undetected 
dwarf host galaxy (M > —12.8). Such a host galaxy 
would be one of the faintest on record; the SN1999aa- 
like SN la SN1999 aw had a host absolute magnitude of 
Mv = -12 .4 (|Strolger et al.ll2002fl , but this was a slow- 
declining luminous event, in contrast to PTF 09dav. Sim- 
ilarly, SN2005E was also located some distance (~ 23kpc 
projected) from the center of its host galaxy, outside the 
plane of the disc, with no evidence for a faint dwarf host 
(or indeed any star formation activity) at the SN posi- 
tion. 

The lack of a detection of a faint host galaxy at the po- 



sition of PTF 09dav argues against a core-collapse origin 
and points to an older progenitor. The likelihood of find- 
ing a high-mass star far out in the galactic halo of a mas - 
sive galaxy is discussed in detail bv lPerets et al.1 (|2010bf ) , 
who conclude that a low-mass, old stellar progenitor is 
substantially more plausible at these large radii. A high- 
velocity and high-mass star ejected from the disk or cen- 
ter of the host of PTF 09dav is always conceivable, al- 
though an 8M Q star would have to be ejected from the 
center at ~ 1000 kms -1 to reach a distance of 41kpc in 
its ~ 40Myr lifetim e. While such stars have been ob- 
serve d in our galaxy (jEdelmann et alJfeOOSl : iHeber et al.l 
2008), they are rare, and represent only a tiny fractio n 
of massive stars in the Milky Way (|Perets et al. 2010b). 
They would also, by themselves, provide no explanation 
for the odd properties of PTF 09dav - the chances of 
PTF detecting a SN from such a star and it simultane- 
ously having peculiar (unique) properties, seem low. Of 
course, the future detection of a faint host for PTF 09dav 
would certainly weight the argument in favor of a core- 
collapsc origin - deeper imaging is underway (Kasliwal 
et al., in prep.). 

At such a large galactic radius, the most likely progen- 
itor type for PTF 09dav is an old (possibly metal-poor) 
white d warf star, perhaps even located in a globular clus- 
ter (e.g. lPfahl et al.ll2009l ). One simple physical explana- 
tion for PTF 09dav is that it is just a fainter example of a 
normal SN1991bg-like SN. The enhanced Ti II compared 
to brighter SN1991bg-like events, together with Sc II and 
Mg I and a lack of S II, point to lower temperatures con- 
sistent with the small inferred Mni, and the light curve 
is well represented by a stretched version of SN1991bg. 
However, the SN also appears quite calcium-rich, partic- 
ularly at later phases, with nebular [Ca II] A7307 begin- 
ning to emerge in the +16day spectrum, and an z-band 
excess consistent with strong [Ca II] emission at day +30 
and beyond. This scenario also provides no immediate 
explanation for strong Sc II in the spectra (other than 
the cool temperatures that would at least allow any Sc II 
present to be seen), and, if the possible He I in the spec- 
trum is real and not just Na, then PTF09dav is a very 
different event to other SN1991bg-like SNe. 

Since calcium-ri ch spectra c o uld aris e from the det- 
onation of He, iPerets et al.1 (|2010bft propose that 
SN2005E-like SNe be explained by the de tonation 
of helium on a CO core w hite dwarf (e.g. iNomot' 
1982 blR IWooslev et all [1981 iWooslev fc Weaverl 1199 
Livne fc Arnettl Il995h . These models are also worth 
exploring in the context of PTF09dav. Various con- 
temporary models of the observable outcomes of these 
sub-Chandrasekhar-mass explosions exist, with a va- 
riety of CO cores, ignition densities, and He shell 
masses modeled, and which consider bot h full star (CO 
core plus shell) and shell-only explosion s (Bildstc n et al.1 
20071 IShen"et al.ll2010t iSim et al.| [2010H Fink et alJ l201C" 



Wald man et al.ll2010tlWooslev fe Kas cn 2010). The 'full 



star' explosions tend to be too luminous to explain the 
properties of sub -luminous SNe, but the helium shell det- 
onation models (jShen et alJl2010b iWaldman et alJl201Ct ) 
are more promising. 

The light curves tend to be powered by the decay of 
56 Ni, 48 Cr, 52 Fe plus other radioactive nuclei and their 
decay products, but the key predicted observables of such 
scenarios can vary significantly according to the configu- 
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ration of the progenitor system. iShen et all (|2010l ) con- 
sider CO cores of > O.6M , and their models are typ- 
ified by faint (Mr = —15 to —18), rapidly rising (<10 
days) and possibly multi-peaked light curves which red- 
den quickly following maximum light due to the emissiv- 
ity of Fe-group isotopes when they transition fr om the 
doubly to singly- ionized state (e.g. Kasen 2006). The 
spectra are dominated by Ca II and Ti II features at 
high velocity (~ 10 4 kms^ 1 ), lacking intermediate mass 
elements (such as silicon). Lo wer CO core masses o f 
< 0.6M© have been explored bv lWaldman et all ([2010), 
who show that at these low masses the abundance of 56 Ni 
drops sharply, with light curves powered predominantly 
by 48 Cr, with the abundances of Ca and Ti both increas- 
ing together with the amount of unburnt helium. The 
presence of unburnt helium in the observed spectra of 
all He detonatio n models is unclea r due to the LTE as- 
sumpt ions in the IShen et all (|2010f ) and IWaldma n et al.l 
(|2010l ) models, but some He might be expected. 

The observations of PTF 09dav match several of these 
features. The spectra contain significant Ca II and Ti II. 
The light curve also peaks with a consistent luminosity, 
and becomes red in (e.g.) r — i at later phases, although 
this may also be plausibly explained by the emergence of 
the nebular [Ca II] doublet at A7307. The ejecta velocity 
and rise-tim e are both on the e dge of the model outcomes 
predicted bv lShen et al.1 (|2010f ) , although more consistent 
with t he lower CO-core mass models of IWaldman et all 
(2010). However, there is no immediate explanation for 
the presence of Sc II (which is not present in any of the 
model spectra) , and the spectra of PTF 09dav also con- 
tain both Si II and O I, which are predicted to have very 
low abundances in these helium detonation events. Thus 
PTF 09dav is unlikely to be the result of a pure helium 
detonation, at least in the forms currently considered in 
the literature. 

Another possibility for PTF09dav is the deflagration 
of a helium shell layer on the surface of a CO white 
dwarf, which may occu r for smaller He shell masses 
(jWooslev fc Kasen|[20T0l) . These can have quite differ- 
ent properties to those of helium shell detonations, with 
the dim light curve (Mb ~ — 15) again powered mostly 
by 48 Cr (rather than 56 Ni) leading to abundant Ti in 
the spectra. As the burning is incomplete, as with the 
detonation models some helium should also remain. The 
explosion energy per unit ejected mass is smaller than for 
the detonation models, giving slower light curves than for 
detonation models of the same brightness, with relatively 
low ejecta velocities ~ 4000 kms -1 . The presence of Sc in 
the spectra (from 45 Ti) is a prediction in t he helium de- 
flagration models ()Wooslev fc Kaser]|2010D . although it 
is dependent on the neutron excess and hence metallic- 
ity. However, these deflagration models also predict quite 
blue optical colors (B — V ~ 0) with peak luminosities 
fainter than PTF 09dav, and, as for the detonation mod- 
els, there is a lack of Si/Mg in the model spectra. In this 
regard, PTF09dav seems to combine both C/O-burning 
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and He-burning products - further detailed modeling and 
comparisons with observations will be required for a full 
understanding. 

Whatever the physical explanation, PTF09dav illus- 
trates the discovery power of new high-cadence transient 
surveys such as PTF. Although the 'faint and fast' obser- 
vational parameter space remains relatively unexplored, 
these early new discoveries are able to test our under- 
standing of cosmic explosions and motivate new direc- 
tions in explosion modeling. 
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Fig. 1. — A false color image constructed from P60 g, r and i data showing the location of PTF09dav in relation to the nearby "host 
galaxy". The upper arrow indicates the putative host, and the lower arrow the SN. The separation is 56.8", which at z cmb = 0.0359 
corresponds to 40.6 kpc. The host galaxy has an apparently disturbed morphology, and a spectrum shows it to be dominated by nebular 
emission lines indicative of ongoing star formation. No host is detected at the position of PTF09dav to mpiG6lO = 23.2 (Mrq610 = —12.8). 
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Photometry for. PTF 09dav. Photometry is listed for a common zeropoint of 30 a . 
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a Convert (positive) counts, c, to Vega magnitudes, m, using m = — 2.51og(c) + 30 
b Modified Julian Date; MJD = JD - 2400000.5 
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TABLE 2 

Photometric properties of PTF09dav compared to 1991bg-like SNe Ia 
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a Corrected for Milky Way extinction but not any host galaxy extinction. Assumes Ho = 
70kms" 1 Mpc" 1 

b Relative to our subluminous SN Ia template. We emphasize that these stretches cannot be 
compared to stretches derived using a normal SN Ia light curve template which have a very 
different definition of s = 1. 

c Estimated at the peak in the bolometric light curve 

d Assumes that T r for each SN is given by the SN stretch in column 4 multiplied by the 
template rise-time of 14 days. 
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Fig. 2. — A comparison of the + 4.9 day W HT spectrum o f PTF 09dav with a selection of othe r sub-luminous events. From top to bottom: 
SN20 02bj (Poznanski et al. 2010), SN2002cx (Li ct al. 2003 ), SN1991 bg ( Filippc nko et alj 19921 ). PTF09dav, SN2005bl fTaubcnberg er et al.l 
12008) . SN2007ax JKasliwal et al.H2008D . and SN2008ha HFolev et aTll2009l ). For each event the spectrum nearest to +4.9 days is shown. All 
spectra have been deredshifted, and additionally the spectra of SN1991bg and SN2005bl have been shifted by — 3000 km s — 1 to account for 
the velocity difference with PTF 09dav. The vertical arrows highlight the principle differences between PTF 09dav and other sub-luminous 
SNe la (see Fig. \3\ for the evolution of these features in later spectra). 
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Fig. 3.— The +11. 2d and +16d spectra of PTF09dav compared to SN2005M at +12.9d. SN2005bl has been shifted by -3000 km s 
account for the velocity difference with PTF 09dav. The vertical arrows mark the same features as in Fig. [2] 
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Fig. 4. — The observed light curve in apparent magnitude space (left) and flux space (right) for PTF09dav in P48 R and P60/FTN g, 
r, i. For clarity the FTN and P60 data are denoted by the same symbol (but the correct filter responses are used in the light curve fit). 
The overlaid lines are the best-fitting SiFTO light curve using a sub-luminous SN la spectral template, with an arbitrary flux-scaling in 
each filter, but a common "stretch" factor and time of maximum light. Only datapoints up to +20 days are used in the fit; the template 
lines are dotted after this phase and show the predicted evolution at later phases according to the template. Maximum light occurs at 
an MJD of 55051.9. The arrows mark the epochs (relative to maximum light in the rest-frame B-band) of the three photospheric spectra 
in Fig. [5] The reduced \ 2 of the fit is 1.02, and in gRr the light curve decay is completely consistent with that of a stretched version of 
SN1991bg-like SN la. Note that at phases later than +30 days, there is excess i flux compared to the simple template prediction. 



16 



Sullivan et al. 



Till 
Sell/Fell 



Sell/Till 
CrII | Mgl/ScII 



ScII/OI Sell/Till Till/OI 
Sell Nal Sill ! Sill ! Till ! Till 



OI Call 




4000 



5000 6000 7000 

Rest Wavelength (Angstroms) 



8000 



9000 



Fig. 5.— The three photospheric epoch spectra of PTF09dav from WHT (+4.9 days), P200 (+11.2 days) and Keck (+16.0 days), along 
with their SYNAPPS fits. We have assumed a heliocentric redshift of 0.0371. The major features in the spectrum are labeled at the minimum 
of their absorption. They include ionic species common in sub-luminous SNe la (O I, Ca II, Si II, Ti II) together with additional lines 
including Sc II, and relatively strong Na I. We also tentatively identify Sr II and Cr II. Over the period in time covered by these spectra, 
the photospheric line velocities evolve from ~ 6100kms —1 to ~ 5100kms _1 to ~ 4600 km s — 1 . 
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Fig. 6. — A comparison of the photometric properties of PTF09dav and other 1991bg-like subluminous SNe la. Top left: stretch versus 
Mg, top right: Mg versus color ((B — V) ma , x ), and lower left: stretch versus (B — V) ma x- Here, the stretch is measured relative to our 
sub-luminous template and is not on the same system as those used in cosmological studies. All magnitudes are corrected for Milky Way 
extinction. In all panels, PTF09dav is shown as the filled circle, and the other subluminous SNe la denoted by the open circles. The data 
are given in Table [2] 
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Fig. 7. — SYNAPPS fits to the +4.9 day WHT spectrum, comparing Sc II and S II (upper left), He I and Na I (upper right), and the 
evidence for Sr II (lower left) and Cr II (lower right). In all panels, the black line is a section of the WHT spectrum, the red line our default 
SYNAPPS fit, and the blue line a fit with the relevant change in the SYNAPPS set-up. Upper left: S II clearly cannot fit features well-fit by Sc 
II, in particular the A5527 and A5658 lines. Upper right: The evidence in favor of Na or He is less clear-cut, although Na provides a more 
consistent fit. Lower left: Sr II improves the spectral fitting around 4000A. Lower right: Cr II improves some regions of the fit (especially 
around 4800A), but degrades the fit in other areas. 



